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It has been recently found that the equations of motion of several semiclassical systems must take 
into account anomalous velocity terms arising from Berry phase contributions. Those terms are for 
instance responsible for the spin Hall effect in semiconductors or the gravitational birefringence of 
photons propagating in a static gravitational field. Intensive ongoing research on this subject seems 
to indicate that actually a broad class of quantum systems might have their dynamics affected by 
Berry phase terms. In this article we review the implication of a new diagonalization method for 
generic matrix valued Hamiltonians based on a formal expansion in power of h. In this approach 
^"V both the diagonal energy operator and dynamical operators which depend on Berry phase terms 

and thus form a noncommutative algebra, can be expanded in power series in h. Focusing on 
the semiclassical approximation, we will see that a large class of quantum systems, ranging from 
relativistic Dirac particles in strong external fields to Bloch electrons in solids have their dynamics 
C ■ radically modified by Berry terms. 

PACS numbers: 

INTRODUCTION 

Since the seminal work of Berry .1] , the notion of Berry phase has found several applications in branches of quantum 
physics such as atomic and molecular physics, optics and gauge theories. Most studies focus on the geometric phase 
that a wave function acquires when a quantum mechanical system has an adiabatic evolution. Yet, the Berry phase 
in momentum space has recently found unexpected applications in several fields. For instance, in spintronics, such 
a term is responsible for a transverse dissipationless spin-current in semiconductors in the presence of electric fields 
0. In optics, it has been recently found that a monopole in momentum space causes the gravitational birefringence 
■ of photons in a static gravitational fields 3]. Both effects are two manifestations of the spin Hall effect which can 
| be interpreted at the semiclassical level as due to the influence of Berry curvatures on the semiclassical equations of 
motion of spinning particles Similarly, a new set of semiclassical equations with a Berry phase correction was 
proposed to account for the semiclassical dynamics of electrons in magnetic Bloch bands [5(p|. In a more exotic 
application, intrinsic Berry phase effects in the particle dynamics of the doubly special relativity theory was recently 
described 0. 

In the above cited examples, the semiclassical equations of motion with Berry phase corrections (anomalous ve- 
locity terms) can be derived in a representation where the Hamiltonian is diagonalized at the semiclassical order. 
It was indeed shown that the semiclassical diagonalization results in an effective energy operator with Berry phase 
corrections as well as noncommutative covariant coordinates and momentum operators. These dynamical operators 
being corrected by Berry terms this leads directly to these new Berry effects. This is another illustration of 
the fact that the physical content of quantum systems is most often best revealed in the representation where the 
Hamiltonian is diagonal. The paradigmatic example is provided by the Foldy-Wouthuysen (FW) representation of 
the Dirac Hamiltonian for relativistic particles interacting with an external electromagnetic field. In this represen- 
tation the positive and negative energy states are separately represented and the non-relativistic Pauli-Hamiltonian 
is obtained 01. Actually even if several exact FW transformations have been found for some definite classes of po- 
tentials 10] [11] 12], the diagonalization of matrix valued Hamiltonian is, in general, a difficult mathematical problem 
requiring some approximations, essentially a perturbation expansion in weak fields. To overcome this limitation we 
have recently proposed a new method based on a formal expansion in powers of the Planck constant h [131 ] which is 
not restricted to Dirac Hamiltonians but also applicable to a large class of quantum systems. It is worth mentioning 
that recently a variant of the FW transformation valid for strong fields and based also on an expansion in h of the 
Dirac Hamiltonian was presented [14j . The main advantage of the diagonalization procedure of [13j is that it embraces 
several different physical systems ranging from Bloch electrons in solid to Dirac particles interacting with any type of 
external fields (for instance in refs. [3j[l5| electrons and photons in a static gravitational field were considered). 
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In this paper we review the recursive diagonalization procedure of ref. [l3j from which can deduce the expressions 
of the semiclassical energy and of the dynamical operators. We then consider, as a physical application, an electron 
in a magnetic Bloch band and the spin hall effect of light in a gravitational field 

RECURSIVE DIAGONALIZATION OF QUANTUM HAMILTONIAN 

In this section we consider a quantum mechanical system whose state space is a tensor product L 2 (7l 3 ) <8 V with 
V some internal space. In other words, the Hamiltonian of this system can be written as a matrix Hq (P, R) of size 
dim V whose elements are operators depending on a couple of canonical variables P and R, the archetype example 
being the Dirac Hamiltonian with V = C 4 . In [13j we found a diagonalization process for this matrix valued quantum 
Hamiltonian Hq (P,R) recursively as a series expansion in powers of H which gives the quantum corrections to the 
diagonalized Hamiltonian with respect to the classical situation h = 0. For example, the first order correction in 
h corresponds to the semiclassical approximation. In this approach we derived the H expansion recursively in the 
following way. The Planck constant h is formally promoted to a dynamical parameter a in order to establish a 
differential equation connecting the two diagonalized Hamiltonians at h = a and h = a + da. The integration of this 
differential equation allows then the recursive determination of the different terms in the expansion of the diagonalized 
Hamiltonian in powers of a. To start with, consider the diagonalization at the scale a 

U a (P, R) Ho (P, R) U+ (P, R) = e a (P, R) if [P, R] = -ia (1) 

and similarly for a + da. 

U a+da (P, R) H (P, R) U+ +da (P, R) = e a+da (P, R) if [P, R] =-*(<* + da) (2) 

Let us develop this last relation to the first order in da, 

e a+da (P, R) = U a H U+ + da (d a U a H Q U+ + U a H Q d a U+) (3) 

After rewriting the r.h.s. of this equation in terms of Berry connections A^ 1 = iU a V p t U+ and A^ 1 = —iUaW^U^ 
we arrive at the following differential equation [l3| : 

d e a (P,R) = [9 Q C/ Q (P,R)t/+(P,R),e a (P,R)] 



do 



| \A* l X7 Rl£a (P, R) + X7 Rl£a (P, R) A* + A%V Pl e a (P, R) + V P ,e Q (P, R) A{ 
\ { [e a (P, R) , Af } A P J - [s a (P, R) , A* ] A* - [e a (P, R) , [A* ,A%]]} 



- {Asym{V Pl V Rl £ a (P,R)} - U a A.sym {V Pi V Rl H (P,R)}C7+} 

~ [Xe Q (P,R)-£ Q (P,R)X+] (4) 



with the notation 



X = (A S ym[V ff V Pl [/ Q (P,R)])C/+(P,R) (5) 
where the linear operation Asym [8| acts on a symmetrical function in P and R in the following way : 

Asym | \A (R) B (P) + \b (P) A (R)} = ~ [B (P) , A (R)] (6) 

the functions A (R) and B (P) being typically monomials in R and P arising in the scries expansions of the physical 
quantities. 

As shown in [l3| , we can separate the energy equation Eq. Q in a diagonal and a non diagonal part such that we 
are led to the following two equations 

-£-e a (P,R) = V+ [R.H.S. of Eq. @] (7) 

da 

= V- [R.H.S. of Eq. U (8) 
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These two equations are supplemented by the differential unitarity condition 

= d a U a (P,K)U+(P,R)+U a (P,R)d a Ui(P,R) -t(X-X+) + ± [A*, A*] (9) 

The three equations Eqs.([7])-(j9]) allow to determine recursively in powers of a the energy of the quantum system in 
question. Actually, the integration over a of Eq. (|7|) gives e a (P,R) at order n in a when knowing all quantities 
at order n — 1. By the same token, Eqs.® and © (whose meaning is that U a (P,R) is unitary at each order in 
a) involve d a U a (P,R), and allow to recover U a (P,R) at order n by integration over a. As a consequence, the 
diagonalization process is perfectly controlled order by order in the series expansion in a. 



THE SEMICLASSICAL APPROXIMATION 



In this section we consider the Hamiltonian diagonalization at the semiclassical level and the resulting equations of 
motion. Actually, the semiclassical approximation has recently found new important applications in particle and solid 
state physics. Notably, the equations of motion reveal a new contribution coming from the Berry curvature. This 
contribution, called the anomalous velocity, modifies profoundly the dynamics of the particles. For instance, the s pin 
Hall effect of electrons and holes in semiconductors [2J, as well as the new discovered optical Hall effect @| 
can be interpreted in this context. Similarly, the recent experimental discovery of the monopole in momentum can 
also be elegantly interpreted as the influence of the Berry curvature on the semiclassical dynamics of Bloch electrons 

MM- 



The semiclassical energy 

The consideration of Eq. (J4| alone is sufficient to deduce the semiclassical diagonal Hamiltonian. Indeed, writing 
£a — £o + ai£i, with £o the diagonalized energy at the zero order, Eq. ^ is solved by (putting a = ft) 

e (P, R) = e (P, R) + h {K ! V R! e (P, R) + V fii e (P, R) Ag ' + <'Vp,£o (P, R) 

+ V Pi £ (P, R) A$ } + f V+ { [e (P, R) , A*'] <" - [e (P, R) , A$] A*' } (10) 

where we have introduced the notations = V+ [A^] and = V+ [Aq] . 
This latter expression can also be written 



e(p,r) ~ £ (p,r) + yP+ 



e (p,r),^ ! 



v4.n 



e (p,r),^ ; 



v4. n 



o(h 2 



where we have defined the projected dynamical operators (covariant coordinates and momentum operators) 



r = R 
p = P 



hAl 



R 



(11) 



(12) 



with Af = i [U V P U+] ,Ao=~i [U VrU+] ,and <' R - [V R VpC/ ] U+. 

The matrix Uq (P,R) is the diagonalization matrix for Hq when the operators are supposed to be commuting 
quantities, the diagonalized energy being £ (P,R) ■ When P and R do not commute, the matrix Uq (P,R) does 
not diagonalize Hq anymore. In order to get the corrections to the energy at the semiclassical order due to the 
noncommutativity of P and R we have to compute £\ (P, R). From the diagonal Hamiltonian, we can now derive the 
equations of motion for the covariant operators. 



The equations of motion 



Given the Hamiltonian derived in the previous subsection, the equations of motion can now be easily derived. The 
evolution equations have to be considered, not for the usual position and momentum, but rather for the projected 
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variables r and p . Actually, these latter naturally appear in our diagonalization process at the h order. Let us remark 
that their components do not commute any more. Actually 



[n,rj 

\pupj 



] = ih 2 eij = ih 2 (v Pi A Rj - v P] a Ri ) + h 2 [A R] ,A Rt ] 



ih 2 QPf 



^ih 2 {V R% A P -V Ri A Pi 



h 2 [A! 



] = -ihSij + ih 2 e% = -ihSij - ih 2 {V Ri A Rj + V Pj A Pi ) + h 2 [A Pi ,A Rj ] 



the &ij being the so called Berry curvatures. 

Using now our Hamiltonian yields directly to general equations of motion for r, p : 

i r / m i \ ih 
^[r, £ (p,r)] + - 



r, y^+[[ £ (P,r),A Rl ]A Pl 



[e {p,r), A Pl ]A Rl ] 



P = t [P,£(P,r)] + % - 



P,y^+ [[e (P^),A Rl ]A Pl - [e (p, r) , A Pl ] A Rl ] 



(13) 



(14) 



The commutators can be computed through the previous commutation rules between r and p. The last term in 
each equation represents a contribution of " magnetization" type and has the advantage to present this general form 
whatever the system initially considered. In the context of Bloch electrons in a magnetic field, it gives exactly the 
magnetization term revealed in [f| (see [f| and below). For spinning particles in static gravitational fields, this term 
gives a coupling between the spin and the intrinsic angular momentum with magneto-torsion fields 0] (lBj . 



PHYSICAL APPLICATIONS 

Electron in a magnetic bloch band 

This topic was first dealt with in [5j in the context of wave packets dynamics. In the context of the Hamiltonian 
diagonalization it was considered in p] and Q . The purpose is to find the semiclassical diagonal energy operator for 
an electron in a periodic potential facing an electromagnetic field. To apply our formalism, consider an electron in 
a crystal lattice perturbated by the presence of an external electromagnetic field. As is usual, we express the total 
magnetic field as the sum of a constant field B and small nonuniform part <$B(R). The Schrodinger equation reads 
(Hq — e<^>(R)) W(R) = £^(R) with Ho the magnetic contribution (<fr being the electric potential) which reads 

H = ( — + eA(R) + e<SA(R)^ + V(R), P = -tfiV (15) 
\2m J 

where A(R) and <5A(R) are the vectors potential of the homogeneous and inhomogeneous magnetic field, respectively, 
and V(R) the periodic potential. The large constant part B is chosen such that the magnetic flux through a unit 
cell is a rational fraction of the flux quantum h/e. The advantage of such a decomposition is that for <5A(R) = 
the magnetic translation operators T(b) = exp(iK.b), with K the generator of translation, are commuting quantities 
allowing to exactly diagonalize the Hamiltonian and to treat <5A(R) as a small perturbation. The state space of the 
Bloch electron in the periodic zone scheme [20( is spanned by the basis vectors of plane waves \n, k) = |k) (g> \n) 
with n corresponding to a band index and k vary in R 3 . The state |n) can be seen as a canonical base vector 
\n) = (0...010...0...) (with 1 at the nth position) such that U + (k) \n) — \u n (k)) with \u n (k)) the periodic part (in 
space) of the magnetic Bloch waves @[ll[. In this representation K |n,k) = k|n, k) and consequently the position 
operator is R =id/<9k, implying the canonical commutation relations [R^,K -] = iSy. 

We first perform the diagonalization of the Hamiltonian in Eq. (|15[) for <5A = by diagonalizing simultaneously 
Hq and the magnetic translation operators T. The diagonalization is performed as follows: start with an arbitrary 
basis of eigenvectors of T. In this basis Ho can be seen as a square matrix with operators entries. Ho is diagonalized 
through a unitary matrix f/(K) which should depend only on K (since U should leave K invariant, i.e., UK.U + = K) 
and whose precise expression is not necessary for the derivation of the equations of motion, such that U HU + = 
£(K) — e<f>(lTR,U + ), where £ (K) is the diagonal energy matrix made of elements £ n (K) with n the band index (i.e. 
the diagonal representation of Hq). 

Now, to add a perturbation 5A(H) as in (@j), that breaks the translational symmetry, we have to replace K in all 
expressions by 

K = K + e M^) (16) 



5 



and as the flux (5B on a plaquette is not a rational multiple of the flux quantum, we cannot diagonalize simultaneously 
its components i*Q since they do not commute anymore. Actually 



h[K\K j ] = -ie£ ijk 5Bk(R) 



(17) 



To do the semiclassical diagonalization we replace U (K) by U ^K^j , so that the non projected Berry connections are 

A.R t ~ zt/V x.U + and Ax t = V u^SAkCtVjA^. From these we can define the nth intraband position and momentum 

operators r„= H+A n and k„ ~ K — eA n (k n ) x 5H(r n )/H + 0(h) with A n — P n (U"VjytJ + ) the projection of the Berry 
connection on the chosen nth Band [tjl . It can be readily seen that the matrix elements of A n can be written A n (k) = 
i (u n (k)| Vk \ u n (k)) (see also ref. [21|for the derivation of the position operator in the diagonal representation). What 
is totally new here is the transformation on the momentum operator k n which get also a Berry connection correction. 

Using our general results of section II, the full Hamiltonian Eq. (fT5f can thus be diagonalized through the trans- 
formation U(K) + 4% [A^, Apt] U(J£) plus a projection on the chosen n-th Band as it is usual in solid state physics 
(the so called one band approximation) and we obtain the energy operator of the n— th band as 



Pn 



U K)HU + K 



= Pn 



S (kj - - [£(K),UV Kz U+] e^ k —^UV K] U 



l -uv Kj u+ [S(K),UV Ki U+] e «**^E) 



S n (k n ) - M(K).SB(r n )+0(H 2 



(18) 



where the energy levels £„ (k n ^J are the same as £ n {K) with k„ replacing K. The magnetization A^(K) 
^"(H ^v^-)> -A(K) x _4(K)) can be written under the usual form [2l| in the (k, n) representation 



j )nn 



("Ak )n'n 



(19) 



We mention that this magnetization (the orbital magnetic moment of Bloch electrons) , has been obtained previously 
in the context of electron wave packets dynamics ■ 

From the expression of the energy Eq. I|18p we can deduce the equations of motion (with the band index n now 
omitted) 



r = dE(k)/hdk~ k x 6(k) 
Kk = -eE-er x SB(r) - MdSB/dr 



(20) 



where [r l ,r J ] = iQ^Qi) with z:, (k) = 9 8 ^4 J '(k) — d^A l (k) the Berry curvature. As explained in Q these equations 
are the same as the one derived in [5( from a completely different formalism. 



Photon in a static gravitational field. 



We now apply our general approach to the case of a photon propagating in an arbitrary static gravitational field, 
where goi = for i = 1,2,3, so that ds 2 — goo{dx ) 2 — g^dx i dx^ — 0. As explained in [3J] the photon description is 
obtained by considering first a Dirac massless particle (massless neutrino) and then by replacing the Pauli matrices 
a by the spin-1 matrices S. Therefore we start with the Dirac Hamiltonian in static gravitational field which can be 
written 

H = ^a.P + ^£ B/j7 rgV + *Jr°% (21) 

with P given by P a —h l a (IL)(Pi + jZgPyF^a 1 ) with h l a the static orthonormal dreibein (a = 1,2,3), T"' 3 the spin 
connection components and e Q /3 7 cr 7 = |(7 Q 7 /3 — 7^7° )• The coordinate operator is again given by R =ihd . Note 
that here we consider the general case where an arbitrary static torsion of space is allowed. It is known [22j that 
for a static gravitational field (which is the case considered here) , the Hamiltonian H is Hermitian. We now want to 
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diagonalize H through a unitary transformation U(P). Because the components of P depend both on operators P and 
R the diagonalization at order h is performed by adapting the method detailed above to block-diagonal Hamiltonians. 
To do so, we first write J? in a symmetrical way in P and R at first order in h. This is easily achieved using the 
Hcrmiticity of the Hamiltonian which yields 

H = ~ (Vs^a.P + P + -a^) + ^-ylfV. (22) 

Using the general expression Eq. we arrive at the following expression for the diagonal positive (we have projected 
on the positive energy subspace) energy representation e : 

ApT KB.a (A fl xp).B 

where we have introduced a field B 1 = ~\P 8 T a ^e aPl with T Q/3<5 = h{ (h la dih k l 3 - h lfi dih ka ) + h la Yf - h lf3 Tf 5 the 
usual torsion for a static metric (where only space indices in the summations give non zero contributions). We have 
also defined in Eq. (f2"3"| 



AI\(r).p^ , Ar,(r).p 



-f^lTj^r^J' (24) 

with the 7-th component of the vector I\ as I\ 7 = e e/ 3 7 r?' 3 (r) and the helicity A = Note that the dynamical 

operators are now 

r = R + ft c 2 ^iL^ (25) 

p = P-ft c 2 (^5)V R P (26) 

Interestingly, this semi-classical Hamiltonian presents formally the same form as the one of a Dirac particle in 
a true external magnetic field [8] 23]. The term B.cr is responsible for the Stern-Gerlach effect, and the operator 
L = (AflXp) is the intrinsic angular momentum of semiclassical particles. The same contribution appears also in 
the context of the semiclassical behavior of Bloch electrons (spinless) in an external magnetic field 6] 5] where it 
corresponds to a magnetization term. Because of this analogy and since T a|35 is directly related to the torsion of 
space through T Q/3<5 = h & k h la h^T^ we call B a magnetotorsion field. 

However, this form for the energy presents the default to involve the spin rather than the helicity. Actually one can 
use the property Ap/2p = ha/2 — (A^xp) to rewrite the energy as 

4 p 2s p 

The semi-classical Hamiltonian Eq. (|27[) contains, in addition to the energy term e, new contributions due to 
the Berry connections. Indeed, Eq. (f2T|) shows that the helicity couples to the gravitational field through the 
magnetotorsion field B which is non-zero for a space with torsion. As a consequence, a hypothetical torsion of space 
may be revealed through the presence of this coupling. Note that, in agrement with [251 ] . this Hamiltonian does not 
contain the spin-gravity coupling term S.Vg o predicted in [24]]. 

From Eqs. (f2"S"|) and (12"6")) we deduce the new (non-canonical) commutations rules 

[p\p>] = ih®% 

[p\ri] = /////••' • /7/(->; ; ; (28) 

where O^ 3 = d^iA v , — d^iAy + [A^A^i] where £, 77. mean either r or p. An explicit computation shows that at leading 
order 



Qij _ fc4 ( S -P)P7 v y a: 

^pp — 2 £ 4 nPa v rjP/3 

(E.p)p 7 , 



0« = hc ^^ VriPahfc*^ (29) 
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From the additional commutation relations between the helicity and the dynamical operators [r, , A] = [p$, A] =0 we 
deduce the semiclassical equations of motion 

f = (1 — Q pr ) V p e + px9 rr 

p = - (1 - Qpr) V r £ + r x O pp (30) 

To complete the dynamical description of the photon notice that at the leading order the helicity A is not changed by 
the unitary transformation which diagonalizes the Hamiltonian so that it can be written A = hp.H/p. After a short 
computation one can check that the helicity is always conserved 



d fhp.T, 



= (31) 



for an arbitrary static gravitational field independently of the existence of a torsion of space. 

Eqs. (|3"0)) are the new semiclassical equations of motion for a photon in a static gravitational field. They describe 
the ray trajectory of light in the first approximation of geometrical optics (GO). (In GO it is common to work with 
dimensionless momentum operator p =fc " 1 k with = uj/c instead of the momentum [l6j]). For zero Berry curvatures 
we obtain the well known zero order approximation of GO and photons follow the null geodesic. The velocity equation 
contains the by now well known anomalous contribution px O rr which is at the origin of the intrinsic spin Hall effect (or 
Magnus effect) of the photon in an isotropic inhomogeneous medium of refractive index n(r) Q [H[ 03 [III • Indeed, 
this term causes an additional displacement of photons of distinct helicity in opposite directions orthogonally to 
the ray. Consequently, we predict gravitational birefringence since photons with distinct helicities follow different 
geodesies. In comparison to the usual velocity r = V p e ~ c, the anomalous velocity term vj_ is obviously small, its 
order v l ± ~ cXV^g^ being proportional to the wave length A. 

The momentum equation presents the dual expression f x <d pp of the anomalous velocity which is a kind of Lorentz 
force which being of order h does not influence the velocity equation at order h. Note that similar equations of motion 
with dual contributions pxO rr and f x O™ were predicted for the semiclassical dynamics of spinless electrons in 
crystals subject to small perturbations [HQ. 



CONCLUSION 



Some recent applications of semiclassical methods to several branches of Physics, such as spintronics or solid state 
physics have shown the relevance of Berry Phases contributions to the dynamics of a system. However, these progresses 
called for a rigorous Hamiltonian treatment that would allow for deriving naturally the role of the Berry phase. In 
this paper we have considered a diagonalization method for a broad class of quantum systems, including the electron 
in a periodic potential and the Dirac Hamiltonian in a gravitational field. Doing so, we have exhibited a general 
pattern for this class of systems implying the role of the Berry phases both for the position and the momentum. In 
such a context, the coordinates and momenta algebra are no longer commutative, and the dynamical equations for 
these variables directly include the influence of Berry phases through the parameters of noncommutativity (Berry 
curvatures) and through an abstract magnetization term. 

Acknowledgment. It is a great pleasure for us to thank Prof. S. Ghosh for collaborations and for having given 
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